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Abstract. A method for determination and two methods for approximation 
of the domain of attraction D a (0) of an asymptotically stable steady state 
of an autonomous, R-analytical, discrete system is presented. The method of 
determination is based on the construction of a Lyapunov function V, whose 
domain of analyticity is D a (0). The first method of approximation uses a 
sequence of Lyapunov functions V p , which converges to the Lyapunov function 
V on D a {0). Each V p defines an estimate N p of D a (0). For any x S D a (0) there 
exists an estimate N p x which contains x. The second method of approximation 
uses a ball B(R) C D a {0) which generates the sequence of estimates M p = 
f~ p (B(R)). For any x S D a {0) there exists an estimate M p x which contains 
x. The cases ||9o/|| < 1 and p(dof) < 1 are treated separately (even though 
the second case includes the first one) because significant differences occur. 
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1. Introduction. Let be the following discrete dynamical system: 

x k+1 =f(x k ) fc = 0,l,2... (1.1) 

where / : — > f2 is an R- analytic function defined on a domain C R™, 6 O and 
/(0) = 0, i.e. x — is a steady state (fixed point) of (|l.lfl . 

For r > 0, denote by B{r) = {x G 1" : ||x|| < r} the ball of radius r. 

The steady state x = of (|l.l(l is "stable" provided that given any ball B(e), 
there is a ball B{5) such that if x G 5(5) then ) fc (a:) G B(e), for fc = 0, 1, 2, ... [T]. 

If in addition there is a ball B(r) such that / fc (cc) — > as fc — > oo for all x G B(r) 
then the steady state x = is "asymptotically stable" £Q. 

The domain of attraction D a (Q) of the asymptotically stable steady state x = 
is the set of initial states x G £1 from which the system converges to the steady 
state itself i.e. 

D a (0) = ft|/ fc (x) "=3? 0} (1.2) 

Theoretical research shows that the D a (0) and its boundary are complicated sets. 
In most cases, they do not admit an explicit elementary representation. The domain 
of attraction of an asymptotically stable steady state of a discrete dynamical system 
is not necessarily connected (which is the case for continuous dynamical systems). 
This fact is shown by the following example. 

Example 1.1. Let be the function / : R — > R defined by f(x) = \x - \x 2 + ^x 3 + 
jx 4 . The domain of attraction of the asymptotically stable steady state x = is 
D a (0) = (-2.79, -2.46) U (-1, 1) which is not connected. □ 

Different procedures arc used for the approximation of the D a (0) with domains 
having a simpler shape. For example, in the case of the Theorem 4.20 pg. 170 
[J the do main which approximates the D a (0) is defined by a Lyapunov function 
V built with the matrix dof of the linearized system in 0, under the assumption 
||(9o/|| < 1- In a Lyapunov function V is presented in the case when the matrix 
dof is a contraction, i.e. ||<?o/|| < 1. The Lyapunov function V is built using the 
whole nonlinear system, not only the matrix dof. V is defined on the whole D a (0), 
and more, the D a (0) is the natural domain of analyticity of V. In this result 
is extended for the more general case when p(dof) < 1 (where p(dof) denotes the 
spectral radius of dof). This last result is the following: 

Theorem 1.1. (see |2jj If the function f satisfies the following conditions: 

/(0)=0 (1.3) 



p(dof) < 1 (1.4) 

then is an asymptotically stable steady state. D a (0) is an open subset of O and 
coincides with the natural domain of analyticity of the unique solution V of the 
iterative first order functional equation 

V{f(x))-V(x) = -\\xf 

V(0) = [ ' 

The function V is positive on D a (0) and V(x) +oo, for any x° G dD a (0) 

(dD a (0) denotes the boundary of D a (0)) or for \\x\\ — > oo. 
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The function V is given by 

oo 

V(x)=J2\\f k ( x )f for any x e D a (0) (1.6) 

k=0 

The Lyapunov function V can be found theoretically using relation l|1.6fl . In the 
followings, we will shortly present the procedure of determination and approxima- 
tion of the domain of attraction using the function V presented in [21 0] • 

The region of convergence Dq of the power series development of V in is a part 
of the domain of attraction D a (0). If Do is strictly contained in D a (Q), then there 
exists a point x° G 8Dq such that the function V is bounded on a neighborhood of 
x° . Let be the power series development of V in x°. The domain of convergence 
D\ of the series centered in x° gives a new part D\ \ (Dq f] D\) of the domain of 
attraction D a (0). At this step, the part D \J Di of D a (0) is obtained. 

If there exists a point x 1 6 8(Dq [J D\) such that the function V is bounded on 
a neighborhood of x , then the domain Dq [J D\ is strictly included in the domain 
of attraction D a (0). In this case, the procedure described above is repeated in x l . 

The procedure cannot be continued in the case when it is found that on the 
boundary of the domain Dq [_] D\ \] ... [J D p obtained at step p, there are no points 
having neighborhoods on which V is bounded. 

This procedure gives an open connected estimate D of the domain of attrac- 
tion D a (0). Note that f~ k (D), k £ N is also an estimate of D a (0), which is not 
necessarily connected. 

The procedure described above is illustrated by the following examples: 

Example 1.2. Let be the / : R — > M defined by f(x) — x 2 . Due to the equality 
f k {x) = x 2k the domain of attraction of the asymptotically stable steady state 

oo j 

x = is D a (0) = (—1, 1). The Lyapunov function is V(x) = ^ x 2 . The domain 

fe=0 

of convergence of the series is Dq = (—1, 1) which coincides with D a (Q). □ 

Example 1.3. Let be the function / : f2 = (— oo, 1) — > fl defined by f(x) — 
e j r (i- e ) x ■ D ue to the equality f k (x) — e k + ^i_ e k^ x the domain of attraction of the 
asymptotically stable steady state x = is D a (0) = (— oo,l). The power series 

oo 

expansion of the Lyapunov function V(x) — \f k ( x )\ 2 m i s 

k=o 

oo oo 

V{x) = J2 ( m ~ !) - ^T' 2 ^ (I- 7 ) 

m=2 k=0 

The radius of convergence of the series (|1.7(l is 



ro = lim 



\ 



(m- l)^e- 2fc (l - e - k ) m - 2 = 1 (1.8) 



fe=0 



therefore the domain of convergence of the series (|1.7fl is Dq = (—1,1) C D a (Q). 
More, V(x) — > oo as x — > 1 and V{— 1) < oo. The radius of convergence of the 
power series expansion of V in — 1 is 



r_i = lim 

m — >oo 



™ e k {e k_ 1)m -2 [{m _ 3)e k + 2] 

\j ^ (2e k -l) m + 2 { ' 
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therefore, the domain of convergence of the power series development of V in —1 is 
D-i = (—2, 0) which gives a new part of D a (0). □ 

Numerical results for more complex examples are given in E] • 

2. Theoretical results when matrix A = dof is a contraction, i.e. \\A\\ < 1. 
The function / can be written as 

f(x)=Ax + g(x) for any x G (2-10) 

where A = dof and g : £1 — > is an R-analytic function such that g(0) = and 

i™ ii-ii ~ u - 

Proposition 2.1. If \\A\\ < 1, then there exists r > such that B(r) C and 
||/(z)|| < /or any x G B(r) \ {0}. 

Proof. Due to the fact that lim ^p^jp = there exists r > such that B(r) C Q 
and 

|| fl (a:)||<(l-p||)||a;|| for any x G B(r) \ {0} (2.11) 
Let be a; G B(r) \ {0}. Inequality 12.1111 provides that 

||/(z)|| = \\Ax + g{x)\\ < \\A\\\\x\\ + \\g(x)\\ < (\\A\\ + 1 - P||)||z|| = || z || (2.12) 
therefore, ||/(x)|| < ||^[[. □ 

Definition 2.1. Let be R > the largest number such that B(R) C Vl and ||/(x)|| < 
\\x\\ for any x £ B(R)\{0}. 

If for any r > we have that B(r) C £1 and \\f{x)\\ < \\x\\ for any x G B(r)\{0}, 
then R = +oo and B(R) = Q = M n . 

Lemma 2.2. a. B(R) is invariant to the flow of system J^l.l)) . 

b. For any x G B(R), the sequence (||/ fc (a;) ||)fceN is decreasing. 

c. For anyp>0 and x e B(R) \ {0}, AV p (x) = V p (f{x)) - V p {x) < 0, where 

v 

V P (x) = J2\\f k (x)\\ 2 for x eCi (2.13) 

Proof. a. If x = 0, then / fc (0) = 0, for any k E N. For x G B(R) \ {0}, we have 
||/(x)|| < ||x||, which implies that f(x) G B(R), i.e. B(R) is invariant to the 
flow of system (I1.1H . 

b. By induction, it results that for x £ B(R) we have f k (x) G B(R) and 
|/ fc+1 (o;)|| < ||/ fe (x)||, which means that the sequence (||/ fc (cc) ||)a:gn is de- 
creasing. 

c. In particular, for p > and x G B(R), we have ||/ p+1 (a;)|| < ||/(^)|| < ||^|| 
and therefore, AV p (x) = ||/ p+1 (a;)|| 2 - ||a;|| 2 < 0. 

□ 

Corollary 2.3. For any p > 0, there exists a maximal domain G p C f2 such that 
G G p and for x G G p \{0}, the (positive definite) function V p verifies AV p (x) < 0. 
In other words, for any p > the function V p defined by l2.1, e Jp is a Lyapunov 
function for on G p . More, B(R) C G p for any p > 0. 

Theorem 2.4. B(R) is an invariant set included in the domain of attraction D a (0). 
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Proof. Let be x G B(R) \ {0}. We have to prove that lim f k {x) = 0. 

k — >oo 

The sequence (f k (x))k^n is bounded: f k (x) belongs to B(R). Let be (/^ (x))j e N 
a convergent subsequence and let be lim f kj (x) = y°. It is clear that y° g B(R). 

J— ><x> 

It can be shown that 

l!/ fc (aOII > Il2/°ll for any k G N (2.14) 

For this, observe first that f kj (x) — » y° and (a;)||)keN is decreasing (Lemma 
I2.2|) . These imply that \\f kj (x) || > ||y°|| for any kj. On the other hand, for any k G 
N, there exists kj G N such that kj > k. Therefore, as the sequence (||/ fe (^)||)fceN 
is decreasing (Xemma l2~2*|) . we obtain that ||/ fe ( a; )ll > ll/^'WII > \\v°\\- 

We show now that y° = 0. Suppose the contrary, i.e. y° ^ 0. 

Inequality Ij2.14|) becomes 

||/ fc O)|| > > foranyfceN (2.15) 

By means of Lemma \'2. 21 we have that ||/(y°)|| < 

Therefore, there exists a neighborhood U f( y o) C B(R) of f(y°) such that for any 
z G UffyO-j we have \\z\\ < \\y°\\. On the other hand, for the neighborhood C//( y o) 
there exists a neighborhood U y o C B{R) of j/° such that for any y G U y o, we have 
/(?/) G Uf( y o). Therefore: 

\\f(y)\\ < || 2^° || for any y G C/,o (2.16) 

As f kj (x) — > y°, there exists j such that f kj (x) G JJ^o, for any j > j. Making 
y = / fe ' (a;) in l(2~Tnjl . it results that 

ll/ %+1 (z)ll = ll/(/ k, (ar))ll < \\y°\\ for j > j (2.17) 
which contradicts (|2.15|) . This means that y° — 0, consequently, every conver- 
gent subsequence of (f k (x))k^n converges to 0. This provides that the sequence 
(/ fe (x))fc e N is convergent to 0, and x G D a (0). 

Therefore, the ball B(R) is contained in the domain of attraction of D a (0). □ 

For p > and c > let be N£ the set 

JV= = { x g n : < c} (2.18) 

If c = +oo, then = fi. 

Theorem 2.5. Let be p > 0. For any c G (0, (p + 1)-R 2 ], the set NE is included in 
the domain of attraction D a (0). 

Proof. Let be c G (0, (p + l)i? 2 ] and x G iV£. Then K p (a:) = f] ||/ fc (^)l| 2 < c < 

fe=0 

(p+ l)i? 2 , therefore, there exists k G {0, 1, ..,p} such that ||/ fc (x)|| 2 < R 2 . It results 
that / fe (a;) G B(R) G A»(0), therefore, x G D o (0). □ 

Remark 2.6. R is obvious that for p > and < d < c" one has NZ C NZ . 
Therefore, for a given p > 0, the largest part of D a (0) which can be found by this 
method is Np p , where c. p = (p + l)i? 2 . In the followings, we will use the notation 
N p instead of Np" . Shortly, N p = {x g fi : V p (x) < (p+l)R 2 } is a part of D a (0). 
Let's note that N — B(R). 

Remark 2.7. If R = +oo (i.e. O = R" and \\f{x)\\ < \\x\\, for any x G K \ {0}), 
then N p = R n for any p > and D a (0) = R n . 
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Theorem 2.8. For the sets (N p ) p ^, the following properties hold: 

a. For any p > 0, one has N p C N p+ i; 

b. For any p > the set N p is invariant to f ; 

c. For any x £ D a (0) there exists p x > s«c/i that x £ N p * . 

p 

Proof, a. Let be p > and x £ N p . Then V p [x) = £ ||/ fe (a;)|| 2 < (p + l)i? 2 , 

fc=0 

therefore, there exists k £ {0,l,..,p} such that ||/ fe (a;)|| 2 < R 2 ■ It results that 
f k (x) £ B(R) and therefore f m {x) £ B(R), for any m > k. For m = p+1 we obtain 
||/f +1 (a;)|| < i?, hence V p+1 (x) = V^x) + ||/ p+1 (a;j |[ 2 < (p+ l)i? 2 + R 2 = (p + 2)R 2 . 
Therefore, x £ N p+ i. 

b. Let be x £ N p . If ||a;|| < R then ||/ m (x)|| < R for any m > (by means 

of Lemma QQ). This implies that = £ ||/ fe (/(a;)) || 2 = E ll/^WII 2 < 

fc=0 fc=l 

(p + l)i? 2 , meaning that f(x) £ N p . 

p 

Let's suppose that ||x|| > R. As x £ N p , we have that V p (x) = ||/ fc (a;)|| 2 < 

fe=0 

(p + 1)R 2 , therefore, there exists k £ {0, 1, ..,p} such that ||/ fe (a;)|| < R. It results 
that f k (x) £ B(R) and therefore f m (x) £ B(R), for any m > k. For m = p + 1 we 
obtain ||/ p+1 (a;) ll < R- This implies that 

V P (f(x)) = V p (x) + \\F +1 (x)\\ 2 - \\xf < (p+l)R 2 + R 2 -R 2 = (p+l)R 2 (2.19) 

therefore f(x) £ N p . 

c. Suppose the contrary, i.e. there exist x £ D a (0) such that for any p > 0, 
x N p . Therefore, V p (x) > (p + 1)R 2 for any p > 0. Passing to the limit for 
p — > 00 in this inequality, provides that V(x) = 00. This means x £ dD a (0) which 
contradicts the fact that x belongs to the open set D a (0). In conclusion, there 
exists p x > such that x £ Np* . □ 

For p > let be M p = f- p {B(R)) = {x £ ft : f p (x) £ B(R)}, obtained by the 
trajectory reversing method. 

Theorem 2.9. The following properties hold: 

a. M p C D a (0) for any p > 0; 

b. For any p > 0, M p is invariant to f; 

c. M p C M p+ i for any p > 0; 

d. For any x £ D a (0) there exists p x > such that x £ M p x . 

Proof, a. As M p = f- p (B(R)) and B(R) C D a (0) (see Theorem EHl) it is clear 
that M p C D a (0). 

b. and c. follow easily by induction, using Lemma 12.21 

d. x £ D a (Q) provides that f p (x) — > as p — > 00. Therefore, there exists p x £ N 
such that f p {x) £ B(R) 1 for any p > p x . This provides that x £ M p for any 
p>p x . □ 

Both sequences of sets (M p ) pg N and (N p ) p ^ are increasing, and are made up of 
estimates of D a (0). From the practical point of view, it is important to know which 
sequence converges more quickly. The next theorem provides that the sequence 
(Mp)pgN converges more quickly than (iVp) pe jj, meaning that for p > 0, the set M p 
is a larger estimate of D a (0) then N p . 

Theorem 2.10. For any p > one has N p C M p . 
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Proof. Let be p > and x G iV p . We have that V p (x) = £ ll/ fc (^)l| 2 < (p + 1)-R 2 , 

fc=0 

therefore, there exists fc G {0,1,.., p} such that ||/ fc (a;)|| < i?. This implies that 
f m {x) G B(R), for any m > fc. For to = p we obtain / p (x) 6 B(R), meaning that 
x G M p . □ 

3. Theoretical results when A = dof is a convergent matrix, i.e. p(A) < 1. 
As p(^4) < ||^4||, the case when A is a convergent matrix 0] is more general then 
the case when A is a contraction treated in the previous section. 

Proposition 3.1. If p(A) < 1, then there exists p £ N* and r p > such that 
B{r p ) C n and ||/ p (x)|| < \\x\\ for any p G + 1, .., 2p - 1} ondie B(rp)\{0}. 

Proof. We have that p(A) < 1 if and only if lim A p — (see 01), which provides 

p — >oc 

that there exists p G N* such that ||^4 P || < 1 for any p > p. Let be p G N* fixed 
with this property. 

The formula of variation of constants for any p gives: 

p-i 

f p {x) = A p x + J2 A p ~ k ~ l 9{f\x)) for all x G fi and p G N* (3.20) 

fc=0 

Due to the fact that for any k G N we have lim ^ 9 ^{i = 0, there exists r%5 > 
such that for any p G + 1, .., 2p — 1} the following inequality holds: 
p-i 

J] ll^-lllsCf^))!! < (1 - Pl)||z|| for x G Bfo) \ {0} (3.21) 

fc=0 

Let be x G B(r p )\{0} and p G {p,p + l, .., 2p-l}. Using (|3~27)|) and (|3~2T|) we have 

p-i 

||/*(x)|| = \\A*x + 1 £ i A*- h - 1 g(f\x))\\< 

k=0 

< \\An\\x\\+^\\A*- k - 1 \\\\g{f k {x))\\< 

k=0 

< (11^11 + 1-11^11)11x11 = 11^11 (3.22) 
Therefore, ||/ p (x)|| < ||x|| for p G {p,p + 1, .., 2p - 1} and x G B(r^) \ {0}. □ 

Definition 3.1. Let be p G N* £/ie smallest number such that \\A P \\ < 1 for any 
p> p (see the proof of Provosition f3H\) . Let be R > the largest number such that 
B{R) C ft and \\f p (x)\\ < \\x\\ for p E {p,p + 1, .., 2p - 1} and x G B(B) \{0}. 

If for any r > we have that B(r) C 51 and ||/ p (x)|| < ||x|| for any p G 
{p,p + 1, .., 2p - 1} and x G B(r) \ {0}, then R = +oo and B(R) = O = R". 

Lemma 3.2. a. for any x G B(R) and p G {p,p + 1, -.,2p — 1} the sequence 
(Wf kp ( x )\\)keN is decreasing. 

b. For anyp>p and x G B{R) \ {0}, ||/ p (x)|| < \\x\\. 

c. For any p > p and x G B(R) \ {0}, AV p (x) = V p (f(x)) - V p (x) < 0, where 
V p is defined by \2.13\) . 
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Proof. a. If x = 0, then f p (0) = 0, for any p > 0. 

Let be x G B{R) \ {0}. We know that ||/ p (x)|| < ||x|| for any p G {p,p + 
1, ..,2p - 1}. It results that f p (x) e B(R) for any p G + 1, ..,2p- 1} . 
This implies that for any k G N* we have ||/ fep (x)|| < and ||/ (fe+1)p (x)|| < 
||/ fep (x)||, meaning that the sequence (||/ fep (x)||)fcgN is decreasing. 

b. Let be x G B{R) \ {0}. Inequality ||/ p (x)|| < ||x|| is true for any p G {p,p + 
l,..,2p-l}. 

Let be p > 2p. There exists q G N* and p 1 G {p,p + 1, .., 2p — 1} such that 
p = qp + p' ■ Using a., we have that f p (x) G B(R) and f qp (y) < \\y\\, for any 
y G B(R), therefore 

\\P(x)\\ = \\f q nf P '(x))\\<\\P'(x)\\< N (3-23) 

c. results directly from 6. 

□ 

Corollary 3.3. For any p > p, there exists a maximal domain G p C f2 suc/i 
i/iai G G p and for any x G G p \ {0}, the (positive definite) function V p verifies 
AV p (x) < 0. In other words, for any p > p the function V p is a Lyapunov function 
for on G p . More, B{R) C G p for any p > p. 

Lemma 3.4. For any k > p there exists qk G N such that 

11^+3)^)11 < 11^(^11 < 11/^(^)11 for any x G B(R) (3.24) 



Proof. Let be k > p. There exists a unique qk G N and a unique pk G {p,p + 
1, .., 2p— 1} such that fc = qkP + Pk- Lemma l3~2*l provides that for any x G B(R) we 
have that f qkp (x) G B(R) and \\f Pk {x)\\ < \\x\\. It results that 

ll/ fe (^)ll = \\f Pk (f qkP (x))\\ < \\f qkP (x)\\ for any x G B(R) (3.25) 

On the other hand, we have (qk+3)p — k+(3p—pk). As (3p—pk) G {p+l,p+2, .., 2p} 
and k > p, Lemma 13 . 21 provides that for any x G B(R) we have that f k {x) G B(R) 
and ||/ 3 P- pfc (x)|| < ||x||. Therefore 

\\f^+^{x)\\ = ||/ 3p - pfc (/ fc W)l! < ||/ fc (as)|| for any x G B(i?) (3.26) 
Combining the two inequalities, we get that 

\\f (qk+3)p (x)\\ < \\f k (x)\\ < \\f qkp (x)\\ for any x G B{R) (3.27) 
which concludes the proof. □ 
Theorem 3.5. B(R) is included in the domain of attraction D a (0). 
Proof. Let be x G B(R) \ {0}. We have to prove that lim f k (x) = 0. 

k — >oo 

The sequence (/ fe (x))fcgN is bounded (see Lemma I3.2JI . Let be (f kj (x))j 6 N a 
convergent subsequence and let be lim f kj (x) = y°. 

We suppose, without loss of generality, that kj > p for any j G N. Lemma 13.41 
provides that for any j G N there exists qj G N such that 

||/ fc+3)p (x)|| < ||/^(x)|| < \\f q ^(x)\\ (3.28) 

As (||/ 9jP (x)||)j C N and (||/^ 9j+3 ^(a;)||)j g N are subsequences of the convergent se- 
quence {\\f qp {x) ||) g eN (decreasing, according to Lemma [3.211 . it results that they 
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are convergent. The double inequality (|3.28j) provides that lim ||/ 9j ^(a;)j| = \\y°\\. 

j— >°° 

Therefore, lim \\f qp (x)\\ = \\y°\\. 
It can be shown that 

\\f k (x)\\>\\y°\\ foranyfc>p (3.29) 
For this, remark that lim ||/ 9P (a;)|| = ||y°|| and (||/ 9P (x)||)„ e N is decreasing (Lemma 

q — >oo 

13.2(1 . which implies that ||/ 9 ^(a:)|| > \\y°\\ for any q £ N. On the other hand, Lemma 
l3~l provides that for any k > p there exists q k such that \\f ( - qk+3 ' )i '(x)\\ < ||/ fc (a;)||. 
Therefore, ||/ fc (a;)|| > \\f {qk+3)f> (x)\\ > \\y% for any k > p. 

We show now that y° = 0. Suppose the contrary, i.e. y° ^ 0. 

Inequality l|3.29|l becomes 

\\f k (x)\\ > \\y°\\ >0 foranyfc>p (3.30) 

By means of Lemma f3. 21 we have that ||/^(y°)|| < \\y°\\- 

There exists a neighborhood Uf p( y o) C B(R) of f*(y°) such that for any z £ 
Ufp^yO) we have < On the other hand, for the neighborhood Upi y o\ 

there exists a neighborhood U y o C B(R) of y° such that for any y £ U y o, we have 
P{y) £ Ufffyoy Therefore: 

\\f(y)\\ < \\y°\\ for any y e t/ y o (3.31) 

As f kj (x) — * y°, there exists j such that f kj (x) £ U y o, for any j > j. Making 
y = f kj (x) in l|3.31|l . it results that 

\\fW(x)\\ = \\f(f k >(x))\\ < \\y°\\ for j > j (3.32) 

which contradicts (|3.3U|I . This means that y° = 0, consequently, every conver- 
gent subsequence of (f k (x))k^n converges to 0. This provides that the sequence 
(f k (x))ken is convergent to 0, and x £ D a (0). 

Therefore, the ball B(R) is contained in the domain of attraction of D a (0). □ 

Theorem 3.6. Let be p > p. For any c £ (0, (p + 1)R 2 ], the set N£ is included in 
the domain of attraction D a (0). 

Proof. Let be c G (0, (p + l)R 2 } and x £ N£. Then V p {x) = J2 \\f h {x)\\ 2 < c < 

(p+l)R 2 , therefore, there exists k £ {0, 1, ..,p} such that ||/ fe (a;)|| 2 < R 2 . It results 
that f k (x) £ B(R) C D a (Q), therefore, x £ D a (0). □ 

Remark 3.7. It is obvious that for p > p and < c' < c" one has N£ C . 

Therefore, for a given p > p, the largest part of D a (0) which can be found by this 
method is Np p , where c p = (p + 1)R 2 . In the followings, we will use the notation 
N p instead of Np p . Shortly, N p = {x £ Q, : V p {x) < {p+l)R 2 } is a part of D a (0). 

Remark 3.8. If R = +oo (i.e. Cl = R n and \\f p (x)\\ < \\x\\, for any p £ {p,p + 
1, ..,2p- 1} and x £R\{0}), then N p = E™ for any p > p and D a (0) = R". 

Theorem 3.9. For any x £ D a (0) there exists p x > p such that x £ N p *. 

Proof. Let be x £ D a (0). Suppose the contrary, i.e. x £ N p for any p > p. 
Therefore, V p (x) > (p + 1)R 2 for any p > p. Passing to the limit when p — > oo in 
this inequality provides that V(x) = oo. This means x £ dD a (0) which contradicts 
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the fact that x belongs to the open set D a (0). In conclusion, there exists p x > 
such that x E Np* . □ 

Remark 3.10. In the case when ||(3o/|| < 1 we have shown that the sequence of 
sets (iVp)pgN * s increasing (see Theorem \2.tftl . 

Open Question: If ||<9o/|| > 1 an d p(dof) < 1; is the sequence of sets (N p ) p ^jq 
increasing? 

For p > let be M p = f- p (B(R)) = {x E fi : f p (x) E B(R)}, obtained by the 
trajectory reversing method. 

Theorem 3.11. For the sets (M p ) p6 n the following properties hold: 

a. M p C D a (0), for any p > p; 

b. M kp C M(fe+i) p for any k E N and p E {p,p+ 1, .., 2p - 1}; 

c. For any a; € i5 a (0) there exists p x > p such that x E M p x . 

Proof, a. As M p = f- p (B(R)) and B(R) C D a (0) (see Theorem ESJ it is clear 
that M p C £> Q (0). 

b. follows easily by induction, using Lemma 13.21 

c. x E D a (0) provides that f p {x) — > as p — > oo. Therefore, there exists p x > p 
such that f p {x) E B(R), for any p > p x . This provides that x E M p for any 
p> p x . □ 

Remark 3.12. The sequence of sets (M p ) pe N is generally not increasing (see Sec- 
tion 4-' Numerical examples, the Van der Pol equation). 

Both sequences of sets (M p ) pe N and (N p ) p£ ^ are made up of estimates of D a (0). 
From the practical point of view, it would be important to know which one of the 
sets Mp or N p is a larger estimate of D a (0) for a fixed p > p. Such result could not 
be established, but the following theorem holds: 

Theorem 3.13. For any p > one has N p c M p+ p. 

p 

Proof. Let be p > and x E N p . We have that V p (x) = £ \\f k {x)\\ 2 <(p+ l)i? 2 , 

therefore, there exists k E {0, 1, ..,£>} such that ||/ fc (a;)|| < R. This implies that 
f k+m (x) E B(R), for any m > p. For m=p-k+pwe obtain f p+p (x) E B(R), 
meaning that x E M p+p . □ 

4. Numerical examples. 

4.1. Example with known domain of attraction. Let be the following discrete 
dynamical system 

x k+1 = ±x k (l + x 2 k + 2y 2 k ) 

y k+1 = ly k (l + xl + 2yl) keN (4 ' 33) 

There exists an infinity of steady states for this system: (0, 0) (asymptotically 
stable) and all the points (x, y) belonging to the ellipsis x 2 + 2y 2 = 1 (all unstable). 
The domain of attraction of (0, 0) is D a (0, 0) = {(x, y) E R 2 : x 2 + 2y 2 < 1}. 

As ||9(o,o)/|| = 5> we compute the largest number R > such that ||/(a;)|| < ||x|| 
for any x E B(R) \ {0}, and we find R = 0.7071. 
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For p = 0,1,2,3,4 we find the N p sets shown in Figure ^ parts of D a (0, 0) 
(N p C iVp+i, for p > Q). In Figure the thick-contoured ellipsis represents the 
boundary of D a (0, 0). 

In Figure |3 the sets M p are represented, for p = 0, 6 (M p C M p +i, for p > 0). 
Note that Mq approximates with a good accuracy the domain of attraction. 

4.2. Discrete predator-prey system. We consider the discrete predator-prey 
system: 

Xk+1 Z T k{1 ~ Xk) ~ XkVk with a = * , b = 1, k G N (4.34) 

J/fc+l — b X kUk 2 

The steady states of this system are: (0,0) (asymptotically stable), (—1,0) and 
(1,-1) (both unstable). 

We have that ||<9(o.o)/ll = \-, and the largest number R > such that ||/(x)|| < 
||x|| for any x G B(R) \ {0} is R = 0.65. 

Figure presents the N p sets for p = 0, 1, 2, 3, 4, 5, parts of -D a (0, 0) (N p C N p+ i, 
for p > 0). The black points in Figure E3 represent the steady states of the system. 

In Figure 0] the sets M p are represented, for p = 0, 6 (M p C M p +i, for p > 0). 
Note that the boundary of Mg approaches very much the fixed points (—1, 0) and 
(1,-1), which suggests that M 6 is a good approximation of D a (0). 

4.3. Discrete Van der Pol system. Let be the following discrete dynamical 
system, obtained from the continuous Van der Pol system: 

Xk+1 = Xk 7 Vk , 2 with a = 2, k E N (4.35) 

2/fe+i = x fc + (1 - a)y fc + ax%y k 

The only steady state of this system is (0, 0) which is asymptotically stable. There 
are many periodic points for this system, the periodic points of order 2, 5 being 
represented in Figure El by the black points. 

We have that ||9(o,o)/|| — 2 but p(9(o,o)/) = 0- First, we observe that for p = 2 
we have that (9(o,o)/)^ = O2, therefore, ||(9( 0i o)/) p || = for any p> p. 

The largest number R > such that ||/ p (a;)|| < ||x|| for p G {p,p+ 1, .., 2p— 1} = 
{2, 3} and x G B{R) \ {0} is 7? = 0.365. 

For p = 2, 3, 4, 5, the connected components which contain (0, 0) of the N p sets 
are shown in Figure El We have that iV 2 C N 3 C JV4 C ^5. 

In Figure the sets M p are represented, for p = 0,6. Note that the inclusion 
M p C M p+ i does not hold. 
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Figure 2. The sets M p , p = 076 for 
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Figure 3. The sets N p , p = 075 for Ull 
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Figure 4. The sets M p , p = 0,6 for (|P1 
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Figure 5. The sets N p , p = 2~5 for 
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Figure 6. The sets M p , p = 07/6 for 
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